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ABSTRACT 

We present new radio measurements of the expansion rate of the Crab nebula’s syn¬ 
chrotron nebula over a ^^SO-yr period. We find a convergence date for the radio syn¬ 
chrotron nebula of CE 1255 ±27. We also re-evaluated the expansion rate of the optical 
line emitting filaments, and we show that the traditional estimates of their convergence 
dates are slightly biased. Using an un-biased Bayesian analysis, we find a convergence 
date for the filaments of CE 1091 ± 34 40 yr earlier than previous estimates). Our 

results show that both the synchrotron nebula and the optical line-emitting filaments 
have been accelerated since the explosion in CE 1054, but that the synchrotron nebula 
has been relatively strongly accelerated, while the optical filaments have been only 
slightly accelerated. The finding that the synchrotron emission expands more rapidly 
than the filaments supports the picture that the latter are the result of the Rayleigh- 
Taylor instability at the interface between the pulsar-wind nebula and the surrounding 
freely-expanding supernova ejecta, and rules out models where the pulsar wind bubble 
is interacting directly with the pre-supernova wind of the Crab’s progenitor. 

Key words: supernova remnants 


1 INTRODUCTION 

The Crab Nebula is one of the most intensely studied ob¬ 
jects in_astro£hysicSj_y^ hj[etains_the powr to surprise us 
(see iBiihler fc BlandfordI l2014l : iHesterl boOSl . for recent re¬ 
views). It is the remnant of a supernova explosion in the 
year CE 1054 , which was witnessed by Chinese and other 
astronomers llStephenson fc Greerjiioosl l . The presently vis¬ 
ible nebula is bright at all observable wavelengths, and con¬ 
tains one of the first known pulsars, PSR B0531-I-21. The 
pulsar’s spin frequency is 30 Hz, and it is slowing down 
at a rate of —3.78 x 10~^° Hz s“^ (iLvne et al.l 12013 1. The 
Crab Nebula is the prototype of a pulsar-powered nebula, 
commonly known as a pulsar wind nebula (PWN), where 
the rotational energy lost by the spinning down of the pul¬ 
sar powers the nebula. The energy input from the pulsar, 
which emerges in the form of a wind of magnetic field and 
relativistic particles, inflates a large bubble of relativistic 
fluid, which emits synchrotron radiation. The synchrotron- 
emitting fluid expands into the supernova ejecta, with which 
it interacts both dynamically and by photoionizing them. 
The synchrotron-emitting bubble and the optical-line emit¬ 
ting filaments of photoionized thermal gas constitute the 
bulk of the presently visible nebula. 

The Crab has recently been discovered to produce sub¬ 


stantial flares at gamma-ray wavelengths, where the emis¬ 
sion at energies > 100 MeV incre ases by more than a factor 
of two on timescales of days (e.g. iBfihler fc Blandfordll2014 . 
and references therein). The origin of these flares is not yet 
well understood. Since they are poorly localized in gamma- 
rays, we obtained radio observations of the Crab following 
a gamma-ray flare in 2012 August to look for a radio coun¬ 
terpar t. Our initial results are reported in iBietenholz et"^ 
(l2015li : we did not find any such radio counterparts. We 
did, however, obtain high-quality radio images of the Crab. 
In the present paper, we use those radio images (from epoch 
2012), and compare them to earlier radio images to more 
accurately determine the expansion rate of the synchrotron 
nebula. 


The PWN inflated by the pulsar continues to expand, 
and the Crab nebula is young and close enough that the ex¬ 
pansion has been directly observed at different wavelengths. 
In fact, not only is the Crab expanding, it has been accel¬ 
erated since the CE 1054 explosion. Several authors have 


measured the p roper motions of the optical filaments over 
the years (e.g., Duncaji 1 19391: W ol1jej [l95^ : lTrimblelll968l : 
IWvckoff fc Murravlll977l : lNugentlll998l i. and from these cal¬ 
culated a convergence date, which is the date on which the 
filaments would have been closest together if their present 
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positions and proper motions are extrapolated backwards in 
time. This convergence date is later than the known explo¬ 
sion date of CE 1054, indicating an accelerated expansion. 

The expan sion rate of the nebula was also determined 
in the radio by iBietenholz et al.l lll99lll . who used VLA ob¬ 
servations taken between 1981 and 1988. This radio mea¬ 
surement gave a determination of the expansion speed of 
the pulsar wind bubble (as opposed to that of the optical 
line-emitting hlaments), and showed that the former was 
expanding at a rate slightly higher than, but within the un¬ 
certainties, of that observed for the filaments. 

In the canonical picture of the Crab, such acceleration is 
in fact expected, because the PWN is expanding into the su¬ 
pernova ejecta, which themselves are still freely expanding. 
Since the ejecta expand with velocity proportional to the 
radius, as the PWN expands into the ejecta, it must move 
ever faster to overtake them, so as long as there is continued 
energy input from the pulsar and the PWN remains within 
the freely-expanding ejecta, the PWN bubble is expected to 
accel erate (see, e.g. [Reynolds fc Chevalierj 1 19841 : IChevalien 
Il984h . 

If SN 1054 was a normal supernova, then it would be 
expected to have released about 10®*^ erg of energy. The 
presently visible Crab nebula, however, contains at most 
10% of this energy. It is mostly thought that the remain¬ 
der of the expected 10®^ erg resides in freely expanding 
ejecta, outside the presently visible nebula. There has as 
yet been no convincing direct detection of this “outer shell” 
of ejecta, despite the fact that it is predict ed to be substan¬ 
tially ionized by the flu x from the PWN llLundavist et akl 
1 198 A IWang et al.ll2013l') . 

Nonetheless, circumstantial evidence does suggest that 
the presently visible nebula is conhned by the unseen mas¬ 
sive ejecta. In this picture, the optically-bright filaments are 
the product of a Rayleigh-Taylor instability between the low- 
density but hig h-pressure syn chrotron-emitting fluid and the 
massive ejecta. iHesteil (l2008l l gives a summary of the argu¬ 
ments for this interpretation. 

We note here, however, that there is an alternate inter¬ 
pretation, which is that SN 1054 was a low-energy electron- 
capture event, with only a small ejected mass, and that 
the presently visible filaments are the result of the PWN 
interacting with the circumstellar mat erial dSmi thi l2013l: 
iTominaea et al.ll2013l : lMoriva et al.ll2014^ . The progenitor in 
this case is expected to have been a super-asymptotic giant 
branch star, in other words one at the upper end of the mass 
range of the asymptotic giant branch, with a mass of around 
8 to 10 althou gh the mass is not well constrained (see 
iMoriva et al.l [20l3 . and references therein). Such progeni¬ 
tors are expected to produce slow, dense winds before the 
SN explosion. 

If the Crab PWN is currently expanding into the freely 
expanding SN ejecta, then it should, as mentioned above, 
be the case that the synchrotron bubble is expanding more 
rapidly than the optic al filaments. The measurements of 
IBietenholz et al] (Il99ll l already suggested that this might 
be the case, but a better measurement of the expansion of 
the synchrotron bubble is required to be sure. Although the 
synchrotron emission is visible in the optical and even in 
the X-ray, such a measurement is most easily done in the 
radio where the synchrotron bubble is most clearly visible. 
We therefore undertook a new and more accurate determi¬ 


nation of the expansion rate of the Crab’s PWN using radio 
images, as well as a re-evaluation of the expansion rate of 
the optical line-emitting filaments determined from proper 
motion measurements. The radio and optical data allow us 
to measure the expansion rate of two different components of 
the nebula, namely the low-density but high-pressure PWN 
and the massive line-emitting hlaments, respectively. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 VLA observations 

For the purposes of measuring the expansion rate of the Crab 
Nebula, we assemble a collection of 5-GHz images taken be¬ 
tween 1981 and 2012. Images made from data taken at a 
specihc time, and thus only using a single VLA array con- 
hguratioiQ would be best-suited for determining the expan¬ 
sion. We concentrate therefore on using images made only 
from B-conhguration data, which was used for the observa¬ 
tions between 1998 and 20lfl However, due to various rea¬ 
sons, for some of the earlier epochs we use images made from 
combinations of several conhgurations, spatially hltered so 
as to approximate the B-conhguration-only images. We will 
note these epochs below. 

Using observations made only in a single VLA array 
conhguration is the equivalent of applying a spatial hlter 
to the images, limiting both the highest spatial frequency 
(highest resolution) and the lowest one (largest recoverable 
size). In particular, for the B-array at 5-GHz, the resolution 
is limited to ~ l'!2 and the largest recoverable structure is 
~30". Although the range of spatial frequencies sampled is 
similar in all our epochs it will not be identical due to slightly 
different observing frequencies, different observing hour an¬ 
gles and the occasional failure of an antenna. In order to 
unify our images as much as possible, we further spatially 
Hlter our images using a Gaussian kernel. First, we use a spa¬ 
tial low-pass filter to reduce the effective resolution to 2'.'0 x 
T.'8 at p.a. 80°, and secondly a high-pass filter with FWHM 
of 20” to exclude poorly sampled large-scale structure. 

As already me ntioned, we use the 201 2 radio obser¬ 
vations described in IBietenholz et al.l |2oT^. We give here 
only a brief description of the observations: there were two 
sessions of VLA observations on 2012 August 20 and 26, 
spaced 6 days apart (observing code 12A-486). We used a 
bandwidth 2048 MHz around a central frequency of 5567 
MHz, with a total of 5 hours per session. The array was in 
the B configuration, resulting in native resolutions of ^1" 
FWHM, but as mentioned we low-pass filtered them to re¬ 
duce the effective resolution to 2('0 x l”8 for our expansion 
measurement. 

We also used data from 2000 Feb. 11 and 2001 Apr. 
17, described in IBietenholz et al] (l2004fl . These observa¬ 
tions were taken with centre frequencies of 4.885 and 4.625 


^ The VLA array configurations changes approximately three 
times per year. 

^ Note that we do not use images deconvolved with a de- 
fault image constructed from earlier observations as was done in 
IBietenholz et akl (|2oT^,[20o3), since that would bias the measured 
expansion rate to that used to construct the default. 
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GHz. Both frequencies were combined for imaging. We fur- 
ther used data f r om 1 998 August 8, which is described in 
iBietenholz et al] ll200lll . and which was taken with the same 
centre frequencies. 


We further used data taken between 1987 May to 1988 
March, again taken at 4.885 and 4.625 GHz and described 


m IBietenholz & Kronberd (I1990IL 

Bietenholz & Kronbera 

(1991 

). Bietenholz et al.l ( 199lf) and 

Bietenholz & Kronbera 

(1992 

). For these epochs, the calibrated visibility data were 


no longer available, and we worked with the deconvolved 
images, which were made with the combined data from the 
B, GD and D configurations using maximum entropy decon¬ 
volution. We again high-pass filter them with a Gaussian 
of FWHM 20". This filtering will remove almost all of the 
information contributed by the D array, and some of that 
contributed by the G-array. We take as the effective date 
of this image therefore a weighted average of those of the 
B and C configurations only, with the B configuration given 
double the weight of the C configuration. Since all three con¬ 
figurations were observed within one year the effective date 
is not overly sensitive to variations in this weighting and our 
results below do not depend on it. 

Finally, we also used data from B-configuration taken 
on 1982 Oct. 16. Since the signal-to-noise ratio of the B- 
configuration data alone was low, we added some data 
from the G configuration, taken on 1981 Nov. 29. Both 
sets of observations were at 4.9 GHz and had a band- 
width of 12. 5 MHz . The se observations a r e des cribed in 
IWilson et al.l lll985h and IBietenholz et al] (1 199 il l. We re¬ 
reduced the archival visibility data, combined the data from 
both configurations, and then made an image using CLEAN 
deconvolution. We apply the same prescription as we used 
for the 1987 epoch to determine the effective date, of tak¬ 
ing a weighted mean of the B and C configuration observing 
dates, with the B-configuration being given twice the weight, 
and again high- and low-pass filtered the resulting image as 
we did the others. 


3 EXPANSION OF THE SYNCHROTRON 
BUBBLE 

We first measure the expansion rate of the synchrotron bub¬ 
ble using the radio images. To illustrate the expansion we 
show a composite of the images from 1987 and 2012 in Fig¬ 
ure [T] Since there are few sharply-defined features from 
which the expansion could be directly measured, we use 
the sa me approach for d e termining the expan s ion as was 
used in IBietenholz! (|2006ll , IBietenholz &: Bartell (l2008l l and 
IBietenholz et al.l (Il991 1. and which was originally developed 
bv lTan fc Gulll ( 1985ll . We repeat a brief description here for 
the convenience of the reader: rather than determining the 
proper motion of individual features, we measure the over¬ 
all expansion by determining the scaling between a pair of 
images by least-squares. Th is is accomplished by using the 
MIRIAD dSault et al.lll99^ task IMDIFF which determines 
how to make one image most closely resemble another, by 
calculating unbiased estimators for the scaling in size, e, the 
scaling and the offset in flux density, A and b respectively, 
and the offsets in RA and deck, x and y respectively, by least 
squares. Our chief interest is in the expansion factor, e, but 
because of uncertainties in flux calibration, absolute posi¬ 


tion, and image zero-point offsets caused by missing short 
spacings and self-calibration and slightly different observing 
frequencies, all five parameters needed to be determined. 

We calculate the expansion from images which were not 
corrected for the primary beam response. Although the at¬ 
tenuation due to the primary beam is appreciable towards 
the edge of the nebula at 5 GHfl the uncorrected images are 
preferred for the expansion calculation because the noise is 
uniform over the image, which is a desirable property for an 
algorithm that computes a least-squares fit over the whole 
image. The use of the uncorrected images should not affect 
our expansion results because the frequency of observation, 
and thus the primary beam response, was almost the same 
at all our epochs and the expansion is small. We re-imaged 
or re re-sampled all the earlier observations to have the same 
pixel spacing (0.27") as the 2012 ones, and also, if required, 
J2000 coordinates. For each pair of images we compute e. To 
reduce biases, we also repeated the calculation but inverted 
the order of the images, which causes 1/e to be determined. 
As our final value of e was the average of the two runs (al¬ 
though in all cases the two values of e were consistent to 
within 0.001). The uncertainty in e is difficult to estimate, 
we conservatively adopt a value of 0.002 (however, we show 
below that the scatter in the derived expansion rate suggests 
a lower uncertainty in e of ~ 0.0009). We give our values of 
e for various pairs of images in Table O We define the frac¬ 
tional expansion rate, flexp, as the percentage increase in size 
of the nebula per year. The weighted mean value of flexp over 
the period of 1982 to 2012, was 0.134T0.005 % yr~^ (at the 
weighted-mean epoch of 1996.6). 

From each pair of images, we can also calculate a “con¬ 
vergence date” for the nebula under the assumption of 
constant-velocity expansion from a single origin. Let ti and 
t 2 be the times of two images measured with respect to the 
time of the explosion, and At = t 2 — ti. For constant ve¬ 
locity expansion, the value of e between ti and t 2 is just 
(ti -I- At)/ti. From any determination of e we can therefore 
calculate the value of ti and thus estimate the explosion 
date. Since the assumption of constant-velocity expansion is 
well known not to hold for the Crab, we term this date the 
“convergence date” or convergence year, which in the Crab’s 
case will be somewhat later than the actual explosion year 
of CE 1054. 

The weighted mean convergence year calculated from 
our measurements in Table [5] is CE 1243.4 ± 9.3, with the 
uncertainty obtained by scaling the input uncertainties to 
obtain xi = 8. We also performed a Bayesian calculation 
to estimate both the convergence year and the true uncer¬ 
tainties in e, which as mentioned above were not reliably 
estimated a priori. We used a uniform prior distribution for 
the convergence year and a Jeffrey’s prior (p(cr) = 1/cr) for 
the uncertainties in e. The posterior distribution of the con¬ 
vergence year was obtained through Markov Chain Monte- 
Carlo (MCMC) sampling. We obtained a consistent value 


® At the approximate radius of the Crab, 180", from the pointing 
centre and 5.0 GHz, the measured brightness 74% of the true value 
because of the falloff in the primary beam response. 

We used the PyMC package, version 2.3, by C. Fonnesbeck, 
A. Patil, D. Huard, and J. Salvatier to implement the Markov 
Chain Monte-Carlo calculation. This package is available at 
http://pymc-devs.github.lo/pymc/README.html 
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Figure 1. A composite image to illustrate the expansion of the Crab synchrotron nebula. On the top left is the image from 2012 Aug. 
26, while on the lower right is the one from 1987 Dec. Both images are 5-GHz VLA radio images, which have been spatially band-limited 
to show only emission on a limited range of angular scales by convolving to limit the resolution to 2'.'0 X I'.'S at p.a. 80° (FWHM), and 
high-pass filtering at 20” (FWHM). Note that the image brightness can have negative as well as positive values because of the high-pass 
filtering. See text for a fuller description of the data and the processing of the images. 


of CE 1240 ± 12, with the estimate of the measurement un¬ 
certainties in e being 0.0009. 

Is there any bias in this determination of the conver¬ 
gence date? We explored the possibility of a bias as follows: 
we made a copy of a radio image, and expanded it by a factor 
of exactly 1.02 (i.e. by 2%), then added noise. For the noise, 
we used Gaussian random noise with the value in each pixel 
being independent, then convolved with the same CLEAN 
beam (2'.'0x I'.'S at p.a. 80°) as our images, so that our added 
noise has same spatial auto-correlation function as that in 
the real images. We then again used IMDIFF to determine 
e between this modified copy and the original image. In the 
absence of any bias, we expect to obtain e = 1.0200. Over 
n = 8 trials we found that the derived average value of e 


was 1.0205 ± 0.00006, suggesting that any bias in e is less 
than 0.00006. We therefore believe that our estimates of e 
are not significantly biased by the presence of noise in the 
images. 

Does the expansion rate of the Crab change with time? 
We performed a weighted fit to our measured values of Hexp, 
taking the time of each rate measurement as being the mid¬ 
point of the two epochs involved. We plot the resulting val¬ 
ues of Hexp in Figure [21 We found that Hexp = +(5 ± 9) x 
10“"^ ti 054 —(0.341±0.004) % yr~^. Note that flexp represents 
the fractional expansion rate, so for constant velocity expan¬ 
sion since CE 1054, diJexp/dt would be —1.1 x 10~^ % yr“^, 
and for constant velocity expansion since the convergence 
date of CE 1250, it would be —1.8 x 10“'* % yr“^. Although 
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Table 1. Images used for the expansion calculation 


Date of observations 

Julian date 

1982 Apr 24==- 

2445084 

1987 Dec 29 

2447159 

1998 Aug 09 

2451035 

2001 Apr 17 

2452016 

2012 Aug 26 

2456165 


^ For this epoch, to increase the signal-to-noise, we combined 
data taken in the B and C configurations on 1982 Oct. 16 and 
1981 Oct. 29, respectively, but used data only at u-v distance 
> 3.3 KA. As the time of these observations, we take a weighted 
mean, with the B-array, which contributes most of the informa¬ 
tion, being given double the weight of the C-array. The resulting 
date is 1982 June 22. 

^ For this epoch, we used an image made from B,CD and D 
configurations, taken on 1987 Nov. 21, 1988 Mar. 14 and 1987 
May 26, respectively. Due to the spatial high-pass filtering which 
will isolate chiefly the information contributed by the B-array 
we use, as we did for the 1982 epoch, a weighted combination of 
the B and C configuration dates, with again the B-configuration 
being given twice the weight. The resulting date is 1987 Oct. 31. 
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Figure 2. The fractional expansion rate of the Crab Nebula 
(fiexp), as calculated from pairs of radio images. Each blue point 
shows a value of the expansion rate derived from one pair of im¬ 
ages, with the horizontal bar showing the time interval over which 
the corresponding rate was calculated, and the vertical one the 


Table 2. Expansion Factors 


First Epoch 

Second Epoch 

Interval 

Expansion Factor 

1982 

1987 

(yr) 

5.52 

1.007 

1982 

2001 

18.82 

1.024 

1982 

2012 

30.18 

1.040 

1987 

1998 

10.61 

1.015 

1987 

2001 

13.30 

1.018 

1987 

2012 

24.66 

1.034 

1998 

2012 

14.05 

1.019 

2001 

2012 

11.36 

1.016 


our results suggest that Hexp is increasing with time, they 
are well compatible with the expected decrease with time 
(to within Icr), and we do not consider the increase with 
time of Hexp significant (see Fig. ID). 


EXPANSION OF THE LINE-EMITTING 
FILAMENTS 


We also re-examined the expansion rate as determined from 
the proper motions of the line-emitting filaments observed 
in the o ptical. We use d the proper motion measurements 
made bv lNugen'O diggsl i , on the basis of four published high- 
resolution optical images of t he Crab, with the first made 
in 1939 and the last in 1992 ||Baadelll942|; [Gingerichlll977l : 
IParkerlFiggsl : IWainscoat fc KormendJ 199^1 ). The details of 


the astrometrical reduction and proper motion measure¬ 
ments are given in iNugenfj (Il998tl . 

Previously, the expansion age of the Crab had been de¬ 
termined by extrapolating the present positions ba ckwards 
using the (presently measured) proper motions (e.g. iNugentI 
Il998l l. and taking the convergence date as the time at which 
the scatter amongst the extrapolated positions was smallest, 
and the corresponding convergence position as the mean po¬ 
sition at that time. This procedure, however, does not prop¬ 
erly take into account the uncertainties in the proper motion 


uncertainty in the derived expansion rate. The green dashed line 
shows the weighted mean of itexp of 0.134 ± 0.015 % yr“^, with 
the shaded area showing the Icr uncertainty on the mean. Near 
the outside edge of the nebula, i?exp= 0.135 % yr”*^ corresponds 
to 2310 km s“^ {8 = 3', D = 2 kpc). 


measurements, and leads to a slightly, but signihcantly, bi¬ 
ased value for the convergence date, as we will illustrate 
below. _ 

We will assume, like others workers (e.g. iDunca nlll939lL 
that the filaments originated at some particular point in 
space, namely the convergence point, at some particular 
time, namely the convergence time, and that they have 
moved at a constant velocity since then. Since we know in 
fact that the supernova was in 1054 CE, a convergence date 
later than 1054 CE implies that the filaments have not in 
fact been moving at constant speed, but rather have been 
accelerated since 1054 CE. A more straightforward analysis 
would be to fix the explosion date at its known value and 
more directly estimate the amount of acceleration. However, 
we chose the more indirect method of estimating the con¬ 
vergence date for easier comparison with earlier results. The 
amount by which the convergence date is later than 1054 CE 
indicates the amount of acceleration that has taken place. 

The bias in the convergence date can easily be illus¬ 
trated with the following example. Imagine that the Nebula 
were static, in other words that the convergence date was 
infinitely far in the past (or the future). The true proper mo¬ 
tions would be 0. The measured values would be randomly 
distributed about zero because of measurement errors. The 
smallest scatter in the extrapolated positions would be near 
the present time, since no matter what the random motions, 
over any length of time they are more likely to move the fila¬ 
ments farther apart than closer together. The date at which 
the extrapolated position scatter is smallest would therefore 
be near the present, not infinitely far in the past where the 
true convergence date is. This bias occurs even if the real 
proper motions are not zero, in the sense that the date of 
smallest extrapolated position scatter is always somewhat 
nearer to the present than the true convergence date. 
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6 Bietenholz & Nugent 


To take proper account of the errors in both the mea¬ 
sured positions and proper motions, we again turn to a 
Bayesian analysis. We retain the hypothesis that the proper 
motions are constant in time, and that the filaments all orig¬ 
inated at the convergence point in space and at the conver¬ 
gence date in time, and it is the convergence point and date 
that we wish to estimate. For any given convergence position 
and date, the present true filament positions and proper mo¬ 
tions are functionally related. We treat the measurements of 
the positions and proper motions as independent, and in our 
Bayesian analysis we also estimate the true present positions 
of the filaments. We take the measured proper motions and 
positions to be Gaussian-distributed about the true values, 
with standard deviations given by the observational uncer¬ 
tainties. 

We take the following prior distributions: uniform with 
a range of CE 500 to 1500 for the convergence year, Gaussian 
with a — 200" centred around the present mean position 
of the filaments for the convergence position, and Gaussian 
with a with a — 10" centred on the measured positions 
for the true filament positions. Note that the measured fil¬ 
ament positions are very unlikely to be in error by 10”, so 
our prior distribution should be minimally informative and 
should therefore have very little effect on the derived poste¬ 
rior distributions. 

We again turn to Markov Ghain Monte-Carlo to esti¬ 
mate the posterior distributions. We obtain a values for the 
convergence year of CE 1091 ± 34, and the convergence po¬ 
sition of RA = 5** 34“ 33=27, deck = 22° 00' 42"52 (J2000). 

These convergence year and position estimates are not 
sensitive to the exact choice of the prior distribution for 
reasonable choices. We emphasize that our estimate of the 
convergence year is obtain ed using exactly the same mea¬ 
surements as lNueeiitl (Il998l l. and the difference between our 
value and the one Nugent obtained (of CE 1130 ± 16) is due 
our Bayesian estimate not being biased towards the present. 
The bias in Nugent’s (and other earlier values of the conver¬ 
gence year) is not large, although we note that the uncer¬ 
tainty in the convergence year from our Bayesian analysis 
is larger than that of Nugent’s estimate, which is because 
our analysis takes the uncertainties in the measured proper 
motions into account. 


5 DISCUSSION 

5.1 Expansion and acceleration of the Crab 

We found above that the fractional expansion rate of the 
Grab’s synchrotron bubble, Hexp, during the period 1982 to 
2012 was 0.134T0.005 % yr“^, (at mean epoch 1996.6). This 
value of flexp im p lies a convergence year of CE 1255 ± 27. 
iBietenholz et al.l (Il99ll ') found flexp for epoch 1987.4 was 
0.133 ± 0.016 and a convergence year of CE 1243 ± 92 for 
the whole nebula, which agrees well with our determina- 
tion (we note that sin ce we used some of the same data as 
IBietenholz et al]ll99ll . our Hexp estimates are not completely 
independent from theirs). 

It is well known that the Grab N ebula ha s been 
accelerated since the e xplosion (e.g., iTrimbk _[l96^ 
Wvckoff fc Murrain 1 19771 : IBietenholz et al.l Il99ll : NugentJ 
1998l i. The explosion epoch of 1054.5 implies that without 


acceleration, j?exp should be of 0.106 % yr“^ at epoch 1996.6. 
Our measured expansion rate for the synchrotron bubble of 
0.134 ± 0.005 % yr“^ suggests therefore an accelerated ex¬ 
pansion and the later convergence date of CE 1255 ± 27. 
If we assume a power-law expansion, with TZpwn oc t™, 
then we can calculate that m = 1.264 ± 0.049. For a spher¬ 
ical pulsar wind nebula expanding into un-shocked super¬ 
nova ejecta, both theory and simulations predict an ap¬ 
proximately p ower-law expans i on, with m in the range of 
1.1 to 1.3 (e.g. Ghevali ^ 1984 jyagjjerSwahraf et al.l[200ll : 


iBucciantini et al.ll2003l : Gaensler fc Slanell2006l ~l. Our results 

are therefore consistent with the theoretical expectations, 
albeit at the higher end of the range of expected values of 


We also determined the convergence year for the optical 
filaments to be CE 1091T34, which implies a powerlaw expo¬ 
nent of m = 1.040±0.039. Since we found m = 1.264±0.049 
for the synchrotron bubble, we can say that the synchrotron 
bubble is experiencing an acceleration which is stronger than 
that of the optic al filam ents by 3.5cr. 

iRudie et al.l ll2008li examined the proper motion of 
knots in the northern jet, at the northern extremity of the 
Grab, and found a convergence date of CE 1055 ± 24. Al¬ 
though this convergence date is earlier than the biased val¬ 
ues obtained by other workers for the filaments in the body 
of the nebula, it is in fact consistent within the combined un¬ 
certainties with the unbi ased convergence da te we obtained 
for the optical filaments. iRudie et al.l (l2008f) conclude that 
the jet, unlike the filaments in the body of the nebula, had 
experienced essentially no acceleration. Since we find a lower 
acceleration for the bulk of the filaments, the difference be¬ 
tween any acceleration experienced by the jet and that ex¬ 
perienced by the remainder of the filaments is less clear, 
and we cannot conclusively say whether or note the jet has 
experienced less acceleration than the body of the nebula. 

The origin of the optical filaments is generally thought 
to be the following: the pulsar outflow blows a synchrotron- 
emitting bubble, whose interior has high pressure but 
low density, into the still freely expanding supernova 
ejecta. The interface between the pulsar bubble and the 
ejecta is Rayleigh - Taylor unsta ble llGhevalier fc Gullll975l : 
iHester et al.lll99^ : lHesteJl2008l i. This picture is supported 
by the fact that the filaments seem to only occur in a thick 
shell ar ound the exterior of th e nebula, but not in the c entral 
region (iLawrence et al.l 1 19951 : iGharlebois et al.l 120101). The 
magnetohydrodynamic simulations of Forth et al. (l2014h 
show that this instability causes “fingers” of ejecta to de¬ 
velop, which then stream downwards into the PWN, while 
“bubbles” of synchrotron-emitting fluid move outward be¬ 
tween them. 

In this model, on average the synchrotron nebula would 
expand more rapidly than the filaments. Our observations 
show exactly this, and therefore give strong support to the 
idea that the presently visible optical filaments are largely 
the result of the Rayleigh-Taylor and other instabilities at 
the interface between the synchrotron-emitting relativistic 
fluid from the pulsar and the massive supernova ejecta. 

In the alternate scenario of SN 1054 b e ing a low-energy 
electron-capture event (e.g. ISmi t3 I 2 OI 3 I : iTominaea et a.1.1 
l2013l : lMoriva et al.ll2014fl . interaction of the supernova shock 
with the circumstellar medium (GSM) was responsible for 
the luminosity of SN 1054 as well as for the presently visible 
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filaments, which contain largely swept-up CSM rather than 
supernova ejecta. The PWN is therefore confined by the 
dense wind from the super-asymptotic giant branch progen¬ 
itor of the supernova. In this scenario one would also expect, 
as we have observed, that the synchrotron nebula expands 
more rapidly than the filaments. However, our measurement 
of the strong acceleration of the synchrotron bubble is at 
odds with this scenario: if this bubble is confined largely 
by a stellar-wind CSM, expected to have speeds of order 
10 km s“^, or less than 1% of the PWN’s present expan¬ 
sion speed, then one would expect that PWN have been 
strongly decelerated since the explosion. The opposite seems 
to be true, as confirmed by our measurement of the relatively 
strong acceleration of the synchrotron bubble. Our measure¬ 
ment of the acceleration of the synchrotron bubble therefore 
rule out any scenario where the PWN is interacting with a 
slowly moving wind. _ _ 

Recently, lYang fc Chevalieil (|201,5l i suggested a modi¬ 
fied version of the canonical scenario above, where the Crab 
was indeed the result of a low-energy 10®° erg) super¬ 
nova. However, although the energy (and mass) of the ejecta 
is lower than those of a supernova of normal energy (10®^ 
erg), the synchrotron bubble is at present still interacting 
with the freely expandin g ejecta rather than with the CSM. 
lYang fc Chevalieri ll2015h find that for a total ejecta mass of 
4.6 Mq (as found by Fesen et al.ll 19971 ). a supernova energy 
of ^10®° erg is required, in which case the synchrotron bub¬ 
ble is still interacting with the inner, flat density-profile, part 
of the freely-expanding ejecta, although it is approaching 
point in the ejecta density profile where the density profile 
becomes steep. 

Our determination of the relatively strong acceleration 
of the synchrotron bubble (near to the canonical value of 
r oc t®^ ®) is consistent with the low-energy supernova sce¬ 
nario of lYang fc Chevalier! (l2015l i. although it does not dis¬ 
tinguish between that scenario and that of a conventional 
(E ~ 10®^ erg) supernova. In the case of a low-energy event, 
the pulsar bubble might be expected, over the next few cen¬ 
turies, to accelerate further once it reaches the steeper por¬ 
tion of the ejecta density profile, and then decelerate as it 
starts to interact with the CSM, which is moving much more 
slowly than the ejecta. 


6 SUMMARY AND CONCLUSIONS 

1. By comparing our radio images from 2012 with earlier 
ones dating as far back as 1981, we conclude that the average 
fractional expansion rate of the Crab nebula’s synchrotron 
bubble over this period is Hexp = 0.135T0.005 % yr~^. This 
corresponds to a convergence date of CE 1226 ± 27, or a 
powerlaw expansion since CE 1054 with i?pwN oc ii-264±o.049 

2. We also re-examined the proper motion measurements 
made for the optical line-emitting filaments. We find a that 
previous estimates of the convergence year were slightly bi¬ 
ased, and we obtain a new, un-biased Bayesian estimate 
CE 1091 ± 34, corresponding to a powerlaw expansion with 
Rfli oc 

3. The synchrotron bubble shows significantly stronger ac¬ 
celeration than the optical line-emitting filaments. This find¬ 
ing give strong support to the idea that the presently visi¬ 
ble optical filaments are largely the result of the Rayleigh- 


Taylor and other instabilities at the interface between the 
synchrotron-emitting relativistic fluid from the pulsar and 
the massive supernova ejecta. 

4. The relatively strong acceleration of the PWN since the 
CE 1054 seems to rule out scenarios where the PWN is con¬ 
fined largely by a slowly moving CSM, but rather requires 
that the PWN be still expanding into the freely expand¬ 
ing ejecta. This in turn requires that the total energy in 
the ejecta is rather larger than the 5 x lO'^® erg in the 
presently visible filaments. It therefore argues against the 
scenario where SN 1054 was an electron capture supernova 
producing < 10®° erg. 
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